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Abstract 

We present a discussion of the precision for the PHOTOS Monte Carlo algorithm, with 
improved implementation of QED interference and multiple-photon radiation. The main 
application of PHOTOS is the generation of QED radiative coiTcctions in decays of any 
resonances, simulated by a "host" Monte Cai^lo generator. By careful comparisons auto- 
mated with the help of the MC-TESTER tool specially tailored for that purpose, we found 
that the precision of the current version of PHOTOS is of 0.1% in the case of Z and W 
decays. In the general case, the precision of PHOTOS was also improved, but this will not 
be quantified here. 
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1 Introduction 



In the analysis of data from high-energy physics experiments, one tries to resolve the ''experiment — 
theory" equation. This no trivial task requires that a lot of different effects be considered simul- 
taneously. From the experimental side, these are mainly detector acceptance and cuts, which 
are dictated by the construction and physical properties of the detector: the shapes of distribu- 
tions may be distorted by, say, misidentification and residual background contamination; these 
effects need to be discriminated in an appropriate and well-controlled way. From the theoretical 
side, all effects of known physics have to be included in predictions as well. Only then can 
experimental data and theoretical predictions be confronted to determine numerical values of 
some coupling constants or effects of new physics (to be discovered). 

A well-defined class of theoretical effects consists of QED radiative corrections. PHOTOS 
is a universal Monte Carlo algorithm that simulates the effects of QED radiative corrections in 
decays of particles and resonances. It is a project with a rather long history: the first version was 
released in 1991 lUl, followed by version 2.0 Q in 1994 (double emission, threshold terms for 
fermions). The package is in wide use ||3l; recently it was applied as a precision simulation tool 
for W mass measurement at the Tevatron ||4| and LEP |l5l|6l|, and for CKM matrix measurements 
in decays of and B resonances (NA48 0, KTeVlHl , Belle ||9|, BaBar EOl and in Fermilab 

El). 

Throughout the years the core algorithm for the generation of 0(a) corrections did not 
change much; however, its precision, applicability to various processes, and numerical stability 
improved significantly. Increased interest in the algorithm expressed by experimental collabo- 
rations (including future LHC experiments) was a motivation to perform a more detailed study 
of the potential and precision of the PHOTOS algorithm; also new functionalities, such as mul- 
tiple photon radiation and better interference corrections, were recently introduced. The main 
purpose of this paper is however not to document the new features of PHOTOS, as this will be 
done elsewhere lfT2lfT3l . but rather to present the results of precision tests of the algorithm. 

Our paper is organized as follows: in section|2lwe present the physical problem and further 
motivation for our studies. In section |3l we recapitulate the history and properties of the PHO- 
TOS algorithm. In sectionlUwe present the testing procedure and necessary tools we employed 
to obtain our main results, which are presented in section |5] Finally, section |3 summarizes 
the paper and gives an outlook to further developments. In Appendix |X| we describe the most 
important features of the new version of the PHOTOS code. 

Some results of the tests and improvements of PHOTOS, relevant to x physics, have already 
been discussed in lfT4l . 

2 Physical problem and developed solution 

Although QED bremsstrahlung in particle decays is one of the most elementary effects in quan- 
tum mechanics, it is not usually considered explicitly; to simplify, calculations are performed 
for inclusive quantities. In fact, it is only the case of a few specific decay channels, where 
exact fixed-order (e.g. 0{qi?)) fully-differential formulae, with spin amplitudes or matrix ele- 
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ments squared, are available in analytical, semi-analytical or Monte-Carlo form. Nonetheless, 
in the analysis of the experimental data, radiative corrections are usually treated together with 
the detector effects (e.g. conversion, detector efficiency) to form the "QED- subtracted" data. 
The control of the uncertainties of such data becomes a weak point: theoretical uncertainties 
appear on both sides of the ''theory = experiment" equation ifTSll . This problem does not seem 
to be so evident for "discovery" experiments or measurements performed on small statistical 
samples: usually the effects of radiative corrections do not exceed a few per cent. An excep- 
tion is radiative corrections for background processes, such as gg ^ tt; t ^ l^ib, where QED 
bremsstrahlung is an essential element for this channel background contribution to Higgs boson 
searches at the LHC in the // — yy channel lfT6lfT7ll . Conversely, for high-precision measure- 
ments (as those performed nowadays), good control of the radiative corrections becomes vital 
not only for the assessment of the overall experimental error of the respective cross-sections or 
branching ratios, but also for the shapes of the distributions. 

The effects of radiative corrections gradually became an important topic in the context of 
such measurements as high-precision measurements of H^-boson properties (see, for instance, 
the combined results of Tevatron runs [4]), or B, D, K meson decays for measurements of 
CKM-matrix coefficients in 5-physics I18J . With increasing statistics of available experimental 
data, QED radiative corrections have become a significant element in the systematic error of the 
measured quantities. 

Strict, systematic calculations performed order-by-order in the perturbation theory are usu- 
ally not the most efficient way of including the effects of bremsstrahlung. To improve the 
convergence of the perturbative expansion, the most popular method in QED is exponentiation, 
a rigorous scheme of reshuffling the dominant terms between orders of expansion. This method 
is useful for the construction of Monte Carlo algorithms as well lfT9ll20ll . In the leading-log ap- 
proximation, partially inclusive formulae exhibit factorization properties of QED, see e.g. ||2T1 . 
A matrix element formula for particle decay accompanied by bremsstrahlung photon emission 
may be factorized to Bom-level terms times the bremsstrahlung factor^ . 

Similarly, a fully differential formula for Lorentz-invariant phase space for particle decay 
accompanied by a number of photons may be expressed in a way that exhibits factorization 
properties. It is used, for example, in the construction of the TAUOLA algorithm for 
II24II . The important property of this parametrization is the full coverage of the phase space and 
exact (i.e. free of approximations) treatment. The price to pay is that the variables describing 
each added photon are defined in individual rest frames, separated by boosts, making the ques- 
tion of the choice of gauge fixing a very subtle point of the theoretical bases of the algorithm 

El. 

It is nonetheless straightforward to take advantage of these factorization properties and ap- 
proximate the fully differential formula for the cross-section in any particle-decay process ac- 
companied by a bremsstrahlung photon, by a product of a cross-section formula that does not 
contain QED radiative corrections times a bremsstrahlung factor. The bremsstrahlung factor 
depends only on the four-momenta of those particles taking part in the decay, and not on the 

'For some cases this factorization property is also present for non-approximated formulae, e.g. for the (9(a) 
ME formula for Z-boson decays as implemented in MUSTRAAL 1221 1231 at fully differential level. In (Q it was 
shown, contrary to the previous expectations, that this factorization has a natural geometrical interpretation as well. 
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actual underlying process! This approximation, which takes into account both real and virtual 
corrections, converges to an exact expression in the soft-photon region of phase space; this fact 
is exploited in the construction of PHOTOS. It exposes well-known infrared and coUinear singu- 
larities. In the final formula, the infrared divergences that originate from expressions describing 
the emission of real and virtual photons may be regularized and cancelled out order-by-order. To 
realize this operation, a technical parameter giving the minimum photon energy is defined; then, 
integration over the directions of photons with energies lower than the cut-off is performed; this 
results in the final formula being free from infrared regulator. Only photons of energies higher 
than the cut-off are explicitly generated^. The PHOTOS algorithm properly treats the coUinear 
region of the phase space as well: the singularities are regulated simply by the masses of the 
charged particles. 

As a result of the studies of the second-order matrix element, in particular for Z /j+/j^yy 
and gg ?fyy (performed already in 1994 lfT6lfT7l ). the iterative properties of the bremsstrahlung- 
adding formula have been identified and a universal photon-emission kernel has been isolated 
and double-photon emission implemented in PHOTOS. However, at that time, the question of 
precision was left out of the discussion. The main application was an estimate of the background 
to the Higgs boson searches in the H ^ yy channel. Only double hard-photon bremsstrahlung 
was of interest, and a precision of about 10% on the cross-section was sufficient. 

In the iterative algorithm for (fixed-number) multiple-photon emission in PHOTOS, the 
probability of photon generation is based (at crude level) on a binomial distribution, i.e. the 
probability for generating and not generating a photon is calculated in each iteration. Encour- 
aged by the precision of PHOTOS with a photon multiplicity higher than 2 (subsequent terms 
in the expansion improved the agreement with exact calculations), we decided to replace the 
binomial distribution with a Poissonian distribution for the number of photons to be generated. 
In this new Poissonian mode, the algorithm is no longer limited to a fixed maximal number of 
photons: the actual number of photons is also generated (at this "crude-generation" level). The 
modification set the algorithm free from the negative-probability problem as well. At a crude 
level, the total probability for single-photon emission in the iterated kernel does not need to be 
smaller than 1, because it is multiplied by the exponent of an equally large, but negative, num- 
ber-^. Because the crude distribution resembles an expansion of an exponent we call this version 
of multiple-photon algorithm an "exponentiated mode" of PHOTOS. Our approach is also close 
to the language of exponentiation as known in QED since Yennie-Frautchi-Suura times ll26l . 
The exponentiated mode has proved to be more stable, allowing us to significantly lower the 
value of the infrared cut-off on the photon energy as well. 

A flexible organisation of iteration and phase-space variables in the PHOTOS algorithm not 
only allows a full phase space coverage, but also seems to introduce some higher-order effects, 
as will be discussed later in the paper (subsection 15.21) . This was achieved by careful studies 
and comparisons with matrix-element calculations ifT^ [TTl without need of any kind of 
phase- space ordering. 

^From a practical point of view this poses no problem: the energetic resolution of calorimeters used in the 
experiments is also limited. 

^This is trivial: if p is the probability in binomial distribution, it must be smaller than 1; however it does not 
need to have in Poissonian distribution given by P„ = ^e^'^p" 
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3 Evolution of the main features of the PHOTOS algorithm 



PHOTOS is an "after-burner" algorithm, which adds bremsstrahlung photons to already exist- 
ing events, filled in by a host generator (which does not take into account the effects of QED 
radiative corrections) and transmitted by means of a standard HEPEVT event record (only the 
information about four-vectors of particles taking part in the process, and the topology of the 
process are needed). PHOTOS adds (with a certain probability) final-state QED radiation in 
"any" decay of particle or resonance, independent of the physics process that was generated. 
As the result of its execution, bremsstrahlung photons are added in a fraction of the events in 
the HEPEVT event record, taking into account event topology and momentum conservation. 

The single-photon version of the PHOTOS algorithm originates from the MUSTRAAL 
Monte Carlo JEl ES, in its part for the Z a'^a'"(y) process. In PHOTOS the full 0(a) 
matrix-element algorithm was simplified so as to isolate the universal kernel responsible for 
photon emission. The original algorithm, although maintaining full 0{a) precision and expos- 
ing factorization properties, was however dependent on the underlying process"*. A downgrade 
was therefore required to make the algorithm process-independent: the interference terms were 
removed from the formula, and later restored in an approximated way for a limited set of pro- 
cesses (decays of neutral particles into two charged particles of the same mass) by introducing 
a Monte Carlo interference weight. The first version of PHOTOS: universal Monte Carlo for 
QED radiative corrections in particle decays was released in 1991 lUJ. The algorithm operated 
on the four- vectors stored, in the HEPEVT event record lESI . by a host generator. 

In 1994, version 2.0 of PHOTOS was released 0. The most important improvement was 
an implementation of the double-photon emission obtained by iteration of the bremsstrahlung 
kernel. This version became widely used in the HEP experimental era Q. In 2000, PHOTOS 
was integrated to the TAUOLA-PHOTOS-F environment to ease code maintaining; the update 
of PHOTOS was released as part of the documentation of that package ll29ll . In 2003, driven 
by the needs of experimental collaborations, a dedicated Monte Carlo weight for W /v(y) 
process was developed ll30l . This new option, implemented in PHOTOS version 2.07, was 
documented in ll3Tl . One of the changes performed at that time, trivial from the physical point 
of view, but turned out recently to be vital, was the change of variables in the calculation of 
interference weight in PHOTOS: from angles to four- vectors. 

Problems with numerical stability may occur when the four-vectors passed in the event 
record have insufficient precision. PHOTOS is very sensitive to rounding errors and momen- 
tum conservation in the event record, especially when the multiple-photon-emission mode of 
operation is used. At each iteration of the photon-emission kernel numerical rounding errors 
are accumulated. The value of the infrared cut-off parameter was raised to protect the algo- 
rithm from being stopped because of numerical instability. It is not the desired method to solve 
the numerical problems of course. Special kinematic-correcting subroutines have been made 
available in PHOTOS since version 2.07, but they are not active in the default setting of input 
parameters. 

In July 2004, we started a systematic study of the PHOTOS algorithm, focusing on its pre- 

^It remains implemented in the first published version of KORALZ 3.8 f27l. with multiple initial-state 
bremsstrahlung and single final-state radiation. 
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cision and possible extension to the multiple-photon emission. It turned out that the algorithm 
stability suffered from a bug; this had no impact on the results, but prevented the photon emis- 
sion kernel to be iterated more that twice. Once the bug had been identified and corrected, the 
algorithm was extended to triple- and quartic-photon emission, then an exponentiated version 
of an iteration routine was implemented. 

In January 2005, driven by needs of 5-physics experiments, a new, universal interference 
weight was implemented, allowing for the calculation of an interference for "any" process (see 
also subsection 15.41) . Although the formula for this interference weight is given already in 
eq. (17) of |l2|, its integration with the PHOTOS code was not trivial until recently; it was 
modification made in 2003 that made the implementation possible for every decay. 

From the point of view of the precision of the obtained results, the PHOTOS algorithm may 
work in three regimes: (1) as a "crude-level" tool for bremsstrahlung generation in decays of 
any particle or resonance, (2) as a precision tool for dedicated decay channels (at present Z and 
W decays), (3) as a precision tool incorporating matrix-element calculations (not exploited so 
far). 

At a "crude level", only the leading-log soft-photon parts of the matrix element are in- 
cluded in the PHOTOS algorithm^. PHOTOS is most often used as a general-purpose tool, 
thus working at the "crude level" of precision. The precision of the results cannot be guaran- 
teed in this working regime; a typical example of use could, however, be the generation of full 
bremsstrahlung phase-space coverage for acceptance studies only. Appropriate comparisons 
with matrix-element calculations or experimental data would have to be performed to deter- 
mine the actual precision. For the particular decay channel, the PHOTOS algorithm could then 
be upgraded to regime (2) or (3), depending on the form of necessary compensating weight. 
At present, the algorithm employs a number of improvements (with the help of the correcting 
weight) and takes into account such effects as the threshold terms for fermions and the impact 
of the spin of emitting particle. Other effects (e.g. the impact of spin of a decaying particle or 
the influence of spin on interference terms) are not at present taken into account. 

The impact of the approximations and missing terms on the precision of PHOTOS, when 
used in regime (2), that is for W and Z decays, is the main subject of this paper. 

4 Test definition 

Our approach to the study of precision of the PHOTOS algorithm presented in this paper is 
based on numerical comparisons of results obtained from PHOTOS with respect to the results 
produced by other reference Monte Carlo generators. The reference generators that we used 
employ formulae based on the full (i.e. non- approximated) fixed-order matrix-element calcula- 
tions with or without exponentiation; their precision level is well established, both by theoretical 
considerations and by comparisons with experimental data. The advantage of such an approach 

^It can be estimated that the physical uncertainty of the PHOTOS algorithm, with double emission and working 
on the decay of particle P into a charged particle ch and the neutral system y,, is then not smaller than ^ or 

(f log ;^)^, whichever is bigger. 
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is obvious: it allows us not only to check the theoretical precision of the approximation used, 
but also verifies that there are no accidental errors in the actual computing code. 

This approach is probably more appealing to the practically-oriented user than to the theoretically- 
oriented one. We leave more profound theoretical studies for the future, to be possibly per- 
formed in the context of feasibility studies of applicability of methods used in PHOTOS to 
QCD. 

We have tested PHOTOS against a limited number of event generators and for a limited 
number of processes. The choice of processes was dictated not only by physic interests, but also 
by the availability of Monte Carlo generators with well-controlled and established precision. 
The following processes were studied: 

• 2° ^ IJ-IJ+: for this process the LEP era generators KORALZ 1201133 and KKMC US) 
could be used. Note that these programs agreed well with practically all experimental data 
of the LEP measurements. The KKMC generator is based on 0{a^) ME calculations 
with spin amplitudes technique and exponentiation. It can also be used in 0{a) ME 
exponentiated mode. In the case of KORALZ, generation at first-order matrix element 
and no exponentiation were available as well. 

• /j+v^(y): for this process the WINHAC [i33J generator is available for multiple 
photon generation in decay. It is based on first-order matrix element and exponentiation. 
Comparison of PHOTOS with a first-order matrix-element generator without exponenti- 
ation can be found in ref . ll30l ; we will not repeat it here. 

• H ^ /j^/j^ (y) : comparison of PHOTOS with a first-order matrix-element generator with- 
out exponentiation can be found in ref. 0?ll : we will not repeat it here. 

• for leptonic x-decays, complete QED first-order generation with the TAUOLA ll24ll gen- 
erator was available for tests; TAUOLA was widely used and compared successfully with 
LEP and CLEO data ll35ll36ll . 

The event generators and physical initializations mentioned above will be referred to as "ref- 
erence generators". They were used to produce reference results used in comparisons with 
PHOTOS. 

The testbed constructed from the event generators discussed above had one very important 
feature: it allowed a test of PHOTOS against several algorithms that differ in physical content 
and precision. Tests started from an exact 0{a) matrix element for single-photon emission 
(non-exponentiated KORALZ and TAUOLA) compared with a single-photon version of PHO- 
TOS), and went on to multiple-photon generators (exponentiated versions of KKMC, KORALZ 
and WINHAC) compared with triple-, quartic- and exponentiated multiple-photon modes of 
PHOTOS. 

To operate, PHOTOS needs events produced by a "host" generator as an input. We exploited 
the possibility of deactivating the bremsstrahlung generation in reference generators, to turn 
them into QED Born-level "host" generators. 

Let us now define the method for automated tests that were used to obtain the results pre- 
sented later in this paper. To facilitate the systematic comparisons, we adopted a specially 
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designed tool: MC-TESTER, described in detail in ref. ll37ll . The principle of an underlying 
test is to analyse the series of events produced by two distinct Monte Carlo event generators (or 
a sequence of several generators combined together) and to extract characteristic distributions 
in an automatic way. In practice, we search for distinct decay channels of the particle of inter- 
est and store histograms with distributions of all invariant masses that can be formed from the 
four-momenta of its decay products. 

For a selected decay process, in cases discussed here, Z, W and x decays, the distributions 
are extracted from the event record in an automated way (thereby limiting manpower for setting 
up the appropriate analysis code, and also the risk of accidental errors) and stored in the output 
files. Two output files (from distinct runs of event generators instrumented with MC-TESTER) 
are then analysed, and the results are presented in a form visualized as a "booklet" of plots and 
tables. The user is given a general information concerning two compared runs of Monte Carlo 
generators, a list of the decay channels with their branching fractions, and the maximum values 
(for each decay channel) of the Shape Difference Parameter (SDP)^. For each decay channel the 
plots of histogrammed values are then included: each plot presents the two distributions from 
two distinct runs, and the curve being a ratio of the normalized distributions; the value of the 
SDP is also printed for each plot. 

The testing approach implemented in MC-TESTER was already very useful in the case of 
validation of TAUOLA package; however, for the purpose discussed here it required further 
adaptation. The problem was in consistent treatment of the arbitrary number of final-state QED 
bremsstrahlung photons that might be present in the event, the ambiguity caused by infrared 
singularities of QED, which is handled differently by the various Monte Carlo programs. 

Our aim was to develop a technique where comparisons would make physical sense, would 
be automatic and independent of the way infrared singularities are regularized in particular 
generators. For our analysis, we defined zero-, one-, and two-photon topologies in the following 
way: we called the event "zero photon" if there was no photon of energy (in a decaying particle 
rest-frame) larger than the parameter ^test- The "one-photon" event had to have one (and only 
one) photon of energy larger than ^test- If there were more than one such photons, we called 
it a "two-photon" event. If there were more than two photons of energy larger than ^test^ we 
considered only the two most energetic ones, and treated the remaining ones as if they had not 
passed the ^test threshold. For all the photons that did not pass the ^test threshold we performed 
an operation inspired by leading-log logic: we added their four- momenta to the momenta of 
outgoing fermions of smaller angular separation. 

We defined two variants for this test definition: testl and test2. The test! was exactly as 
explained above. In testl, only one photon (the most energetic one) could be accepted^. The 

^The Shape Difference Parameter defined in 1371 quantifies the difference in shape of the histograms coming 
from two compared runs. The histograms contain distributions of all possible invariant masses, which can be 
constructed from the momenta of the decay products of the particle under study; hence there is a set of histograms 
for every decay channel. The SDP value is calculated for each histogrammed mass: it quantifies the exclusive 
surface between the (normalized to unity) corresponding histograms obtained from the two runs. The effects of 
statistical fluctuations are appropriately subtracted. The maximum of SDP over all distributions for a given decay 
channel is taken. 

^If needed, this scheme could be generalized, for instance, to define testS, where up to three photons of energies 
above istest could be accepted, leading to a fourth distinct decay channel. 
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free parameter of the test: Ei^st was adjusted for each process so that the results had physical 
sense. 

Because of the space limitation we cannot present complete "booklets" with results of com- 
parisons; for the purpose of this paper we decided to extract the most vital information: the 
branching ratios for events with 0, 1 and 2 photons (0 and 1 if testl was used) and the maximum 
value of SDP throughout all the plots in a channel. Complete results with summary tables and 
all analysis booklets are available on the web ll38ll . The tests were based on high-statistics runs 
(10^ non-weighted events) for all generators. 

Finally, let us stress that an essential preliminary step for the tests presented here was to as- 
sure the numerical stability of PHOTOS. As a consequence, we could determine that PHOTOS 
may be used for processes at very high energies (i.e. at the scale of LHC energies and even at 
astrophysics ones). The infrared cut on the photon energy may be lowered to 10^^ from the 
10^ value that had to be used before. 

5 Numerical results 

The results presented here cover two issues: precision of the predictions obtained for the single- 
photon emission kernel and the convergence of the iterative solution. At first, in subsection l5.lL 
we will discuss the comparisons of a single-photon emission kernel implemented in PHOTOS 
with other genuine first-order generators. In subsection l5.21 we will discuss tests for multiple- 
photon generation. Finally, in subsection l5.3l we will discuss the applicability of PHOTOS at 
very high energies. 

As mentioned before, we aim at comparing, in total, six versions of distinct Monte Carlo 
programs with five distinct modes of executing PHOTOS. Let us specify them here (including 
acronyms that will be used throughout the text): 

• KORALZ 0(1): KORALZ generator ^ with 0{a) ME for 2° A'^A'"(y), no expo- 
nentiation, single-photon emission 

• KORALZ : KORALZ generator El) with ^(a^) ME for 2° /J^/J^O) and exponenti- 
ation (multiple-photon emission) 

• KKMC : KKMC generator fT^ with 0{a^) ME for P'^P'ij) and exponentiation 
at spin-amplitudes level 

• KKMC 0(1)exp : KKMC generator CD with 0(a) ME for 2° ^ A'^A'"(y) and expo- 
nentiation 

• TAUOLA : TAUOLA generator lEl with 0(a) ME for x /V/Vx(y) (single-photon 
emission) 

• WINHAC EXP : WINHAC generator 01 with full 0(a) ME for W /V/(y) and expo- 
nentiation 

• PHOTOS 0(1): PHOTOS, no iteration 
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• PHOTOS 0(2): PHOTOS algorithm iteration up to two times 

• PHOTOS 0(3): PHOTOS algorithm iteration up to three times 

• PHOTOS 0(4): PHOTOS algorithm iteration up to four times 

• PHOTOS EXP : PHOTOS algorithm with exponentiation, multiple-photon generation 

Before we proceed with a presentation of the results, let us, however, briefly describe the 
approach taken for this presentation and their numerical quantification: 

• As described in section IH we quantify the difference between results produced by two 
generators by calculating branching ratios for events with 0, 1 or 2 photons with energy 
above the Etest threshold, and the maximum values of the SDP parameter of all combi- 
nations of invariant masses in the specific channel (with 0, 1 or 2 photons) for a given 
comparison. 

• In discussion of the results we will also often use a single value of the "overall difference", 
being either the maximum of differences in branching ratios or the maximum of products 
of SDP and corresponding branching ratios^. We then take the larger of the two. 

Later in this section we present only the simplified summary tables and plots where the 
differences in distributions (quantified by the SDP parameter) are most significant. In the sum- 
mary tables the first line (appearing in bold font) refers to the results produced by the reference 
generator, the others refer to results of runs of generators being tested. 

5.1 Single-photon emission kernel 

In this subsection we focus on the quality of the single-photon emission kernel in the PHOTOS 
algorithm. We perform a systematic comparison of PHOTOS with generators that implement 
0(a) matrix element for Z, W and x decays. 

The purpose of these tests was to assess the impact of the simplifications of the photon- 
emission kernel implemented in PHOTOS on the precision of its predictions (in the PHOTOS 
kernel the non-leading 0(a) terms are omitted). 

At first, let us present the comparisons made in the Z A'A'(y) channel. In fact, the original 
PHOTOS algorithm was created by gradual simplification of the exact ME formula for this 
process, so the universal emission kernel could have been identified. Therefore this process 
remains a benchmark for PHOTOS single-emission algorithm, its precision and design. The 
approximations are well understood and may in fact be restored for this particular case. 

We compare PHOTOS with KORALZ Monte Carlo in single-photon emission mode. This 
test, which was already performed a long time ago (when PHOTOS was being developed), has 
been reproduced now, and confirmed that no accidental error has occurred since then. The re- 
sults of this comparison are presented in Fig.[ll One may notice that the difference in branching 

**For instance, in case of the table in Fig.f^the difference in branching ratios is 0.82514 — 0.82362 = 0.00152, 
the product SDP x BR = 0.0053 * 0.176 = 0.00093, so that the first value is taken as the overall difference for that 
test: 0.152%. 
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Figure 1 : Predictions of KORALZ (with 0(a) matrix-element and single-photon emission) are 
compared with predictions of PHOTOS (running in single-photon option) for the l-t^l-t (y) 
channel, ^test =1-0 GeV. The plot presents the distribution of invariant mass of fj fj^ pair 
coming from KORALZ (in red, or darker-grey) and PHOTOS (in green, or lighter-grey); the 
black line is the ratio of the two normalized distributions. The red and green lines are hard to 
separate - they practically overlap. 



GENERATOR 
n photons 


Branchi 



ng ratio 
1 


Max SDP 
1 1 


KORALZ 0(1) 


0.82514 


0.17486 




PHOTOS 0(1) 


0.82362 


0.17638 





0.0053 



Comparison of Mass(1) of mu- mu+ in channel ZO => mu- mu+ gamma | gQp 

0.00534 

1 — 



0.8 



0.6 




MO 
5000 
4000 
3000 
2000 
1000 



20 40 60 80 100 120 



ratios for channel with and without photon are at the per mil level. The differences in distri- 
butions (the maximum value of SDP multiplied by branching ratios) indicate agreement better 
than per-mil: we have quantified the overall difference as 0.15%. 

The dominant contribution to the difference originates from a very sparsely populated area 
of the phase space, where the invariant mass of the lepton system is small, and very hard pho- 
tons were emitted. The discrepancy between the distributions does not exceed 20%, even in 
this comer of the phase space, though. It is visible in the slope of the black curve denoting 
the ratio of the two distributions. In this region of the phase space the approximations (with 
respect to 0(a)) used in PHOTOS are the largest. The differences are too small to justify the 
implementation of a channel-dependent correction weight. 

For the case of -boson decays, we address the reader to QUI , where the comparison of 
PHOTOS and a 0(a) ME generator SANC JSl'^ for the W iuV,nu(y) process is given. This 
paper provides also the theoretical background for the W interference weight implemented in 

ANC calculates complete one-loop amplitudes for the decays of on-shell vector bosons: W, Z or the Standard 
Model Higgs boson; an exact, single, real-photon emission matrix element is obtained that way. 
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Figure 2: Predictions for exact 0{a) in TAUOLA are compared with single-photon emission of 
PHOTOS in the A'~v^Vx(y) channel, ^test = 0.05 GeV. The plot presents the distribution of 
invariant mass of the v^VxY system coming from TAUOLA (in red, or darker-grey) and PHOTOS 
(in green, or lighter- grey); the black curve is the ratio of the two normalized distributions. 
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2003 in PHOTOS 2.07. Presented tests indicate very good agreement between PHOTOS and 
SANC: up to a level of 1% (statistical error) in the majority of the areas of the phase space, 
within 5% for parts of distributions where coUinear-induced logarithms are absent and 10% in 
the regions where only non-leading corrections contribute to the matrix element. 

The paper presents the comparison of various distributions from PHOTOS and SANC 
generators for Higgs boson decays: an agreement at the level of 1% (statistical error) was found 
all over the phase space. 

We consider the results of these tests, published in relatively recent papers ll30l l34l . as 
being relevant and complementary to the tests presented here. Therefore we do not repeat the 
comparison of WINHAC generator ll33l and PHOTOS in single-photon emission mode. 

We have however exploited the opportunity of having another comparison with full 0{a) 
Matrix Element generator, that is TAUOLA ll24l for leptonic x decays. The results of com- 
parison in x^ f^ '^/j'^iiy) decay channel are presented in Fig. |2| The overall difference was 
quantified as 0.11%. The agreement between PHOTOS and TAUOLA results is excellent. In 
the x^ e^VeVx(Y) channel, the agreement is equally good. 
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Table 1: Summary of multiple-photon comparisons in n{y) channel {7r at reso- 

nance peak, CMS rest-frame); iStest = 1-0 GeV. KKMC was used as a reference generator. For 
KKMC-KORALZ comparison the overall difference is 0.00066, for KKMC-PHOTOS EXP it 
is 0.00081, for KKMC-KKMC 0(1)exp it is 0.00137 (all differences are for testl) 
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Max SDP 
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test2 
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test2 
1 
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KKMC 


.83918 
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.01266 




KORAL Z 


.83984 


.16016 


.83984 


.14771 


.01244 
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.0012 


.0012 


PHOTOS 0(2) 


.83925 


.16075 


.83925 


.14630 


.01445 





.0067 





.0035 


.0122 


PHOTOS 0(3) 


.83832 


.16168 


.83832 


.14889 


.01280 





.0038 





.0025 


.0080 


PHOTOS 0(4) 


.83836 


.16164 


.83836 


.14871 


.01293 





.0040 





.0027 


.0058 


PHOTOS EXP 


.83837 


.16163 


.83837 


.14868 


.01295 





.0041 





.0023 


.0092 


KKMC 

0(1)exp 






.83781 


.14881 


.01338 









.0099 


.0467 



5.2 Iterating emission kernels 

In this subsection we cover the comparisons of PHOTOS with other generators for multiple- 
photon generation. In multiple-photon mode, PHOTOS algorithm iterates the single-photon 
emission kernel, which precision was established in the previous subsection. The main ques- 
tions to be answered by the group of tests presented here were the ones concerning the con- 
vergence of the PHOTOS iterative algorithm for photon emission to solutions implemented in 
other generators. 

The KKMC Monte Carlo program [19| played an important role in multiple-photon com- 
parisons: it is the only available Monte Carlo that implements a complete 0{<J?) matrix element 
for two hard photon emission in the Z n^n^n{y) channel. Because of its superior quality with 
respect to KORALZ ll32l (which implements 0{qi}) ME with more approximations) it was used 
as the reference generator throughout all tests in this channel. Before continuing with tests of 
PHOTOS let us first compare the KORALZ and KKMC Monte Carlo programs to assess the 
impact of the presence of second-order terms in the calculations. The results of comparisons 
performed using testl are presented in Fig. 13 The overall difference is below a per-mil level; 
however, with samples of 10^ events the contribution from those O(a^) terms that are missing 
in KORALZ is already visible. As will be shown later, the differences are in fact at the level of 
(and are comparable in shapes to) the differences observed between KORALZ and PHOTOS. 

We start the tests of PHOTOS from comparisons for the benchmark decay Z ij^ij^n{y). 
The complete summary of results is presented in Table[l](as usual, complete results, with MC- 
TESTER booklets are available from the web page ll?Sll ). We tested PHOTOS running with 
various options for photon multiplicity: fixed (two to four) order and the exponentiated versions 
against the reference KKMC generator. 

We observed that adding subsequent iterations of photon-emission kernel improves the 
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Figure 3: Comparison of predictions from KKMC (O(a^) ME with exponentiation at the 
level of spin amplitudes) and KORALZ (O(a^) ME with exponentiation) in the TP /j+/j^(y) 
channel, Etest = 1-0 GeV, testl used for analysis. The plot presents the distribution of invariant 
mass of /j^fj^ pair coming from KKMC (in red, or darker- grey) and KORALZ (in green, or 
lighter- grey); the black line is the ratio of the two normalised distributions. The effects due 
to different types of exponentiation are small, albeit noticeable, and constitute 0.066% overall 
difference. 
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agreement between PHOTOS and KKMC. The difference in results between quartic-iteration 
and exponentiated version of PHOTOS are negligible; however, the exponentiated version is 
technically superior, because it may work with much lower value of the infrared cut-off param- 
eter (see the appendix for details). 

It is striking that the agreement between the exponentiated versions of PHOTOS and KKMC 
is best if the full 0{a^) exponentiated matrix element is used in KKMC; if KORALZ with 
exponentiated 0{a) is used^^ or the matrix element in KKMC is downgraded to exponentiated 
0(a) only (see the last row in the table), the agreement is not as perfect, but it is still good 
enough for any application we can imagine at present. 

The dominant differences are visible in the branching fractions for channels with distinct 
numbers of photons with energies above J^test- The pattern of differences may indicate that 
leading terms summed by an iterated PHOTOS kernel give better precision than any first-order 
exponentiated ME generator for this process. 

To investigate this hypothesis further, we have turned our attention to acoplanarity distribu- 
tions obtained from PHOTOS and KKMC. One of the effects brought by NNLO terms of <9(a^) 
ME is the asymmetry in the acoplanarity distributions. To define acoplanarity, the two planes are 
spanned on momenta vectors of /j and two hardest photons. In Fig.lHwe show acoplanarity 
distributions from the KKMC generator with 0(0^) and 0{a) exponentiated matrix-element 
modes^^. As expected, the distribution is fiat for the 0(a) exponentiated mode and asymmetry 
appears for the 0{a^) exponentiated mode only. 

In Fig. 121 we present the acoplanarity distributions from the exponentiated 0{a^) ME mode 
of KKMC and exponentiated version of PHOTOS. Surprisingly, PHOTOS seems to reproduce 
the bulk of this NNLO effect! This rather subtle effect requires investigation. PHOTOS does 
not use any kind of phase-space ordering, but, at the time when the second photon is generated, 
the momentum of the previous one is changed and so the correlations between the directions of 
two photons appear. We shall not elaborate further on this effect here; however, we believe that 
further enhancements of precision of the PHOTOS algorithm in this aspect may be possible, 
for instance by introducing an asymmetry in the generation of a photon polar angle (uniform 
distribution is currently used in emission kernels). 

Let us now present the results of comparisons of exponentiated mode of PHOTOS and full 
0(a) exponentiated ME predictions of WINHAC for W — > ^V/(y) channels. Both programs 
were running in exponentiated mode. The complete summary of results is presented in Table 
121 The details of one of these comparisons are presented in Fig. |6l The overall agreement is 
better than per mil. The maximum value of SDP in configurations with two hard photons are 
larger for some tests (reaching the level of 0.0255 in the last row of the table); however, this is 
a channel with rather small branching fraction. 

Finally, let us turn to /j^VxV^(y) decay. The comparison of predictions from the 

exponentiated version of PHOTOS and full 0(a) ME of TAUOLA are presented in Fig.|7l The 
overall agreement is better than 0.02%. 

It is interesting to compare the results presented in Fig. |71 with the results in Fig.|2l where 

'"These results are not included in this paper 

"in order to make the effect visible, we used the following cuts: only the events with both photons in the same 
hemisphere as /j^ and having pt > 1.5 GeV were recorded in the histogram. 
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Figure 4: Acoplanarity distributions (in the /j^/j+yy channel) produced by the KKMC 

generator running in exponentiated 0{a^) mode (in red, or darker-grey) and exponentiated 0(a) 
mode (in green, or lighter grey); the ratio of the two distributions plotted in black. For more 
details see subsection l5.2l 




Figure 5: Acoplanarity distributions (in the l-i l-i^YY channel) produced by the KKMC 
generator running in exponentiated 0{a^) mode (in red, or darker-grey) and the exponentiated 
PHOTOS algorithm (in green, or lighter grey); the ratio of the two distributions is plotted as 
dotted black curve. For more details see subsection l5.2[ 
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Figure 6: Multiple-photon comparison of WINHAC and PHOTOS (both running in the expo- 
nentiated mode) in the ^'^\^n{y) channel, ^test = 1-0 GeV, analysed using testl. The 
plot (of largest SDP) presents the distribution of invariant mass of the /j+y pair coming from 
WINHAC (in red, or darker-grey) and PHOTOS (in green, or lighter-grey); the red and green 
lines practically overlap. The black line is the ratio of the two normalized distributions. 
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Table 2: Summary of comparisons of the exponentiated algorithm of PHOTOS and WINHAC 
0{a) exponentiated for leptonic W decays. For more details, see the text. 



GENERATOR 


Branching ratio 


Max SDP 


n photons —>■ 


testl 
1 1 





test2 
1 


2 





testl 
1 





test2 
1 


2 


W+ A'+v^(y) at W mass, £test = 1-0 GeV 


WINHAC EXP 


.92771 


.07229 


.92771 


.07007 


.00222 




PHOTOS EXP 


.92748 


.07252 


.92748 


.07016 


.00236 





.0029 





.0025 


.0023 


W+ A'+v^(y) at W mass, ^test = 5.0 GeV 


WINHAC EXP 


.96491 


.03509 


.96491 


.03473 


.00036 




PHOTOS EXP 


.96470 


.03530 


.96470 


.04488 


.00042 





.0060 





.006 


.047 


W+ e+Veiy) at W mass, ^test = 1-0 GeV 


WINHAC EXP 


.86196 


.13804 


.86196 


.12943 


.00862 




PHOTOS EXP 


.86205 


.13795 


.86205 


.12909 


.00886 





.0022 





.0019 


.0094 


W+ e+Ve{y) at W mass, ^test = 5.0 GeV 


WINHAC EXP 


.93083 


.06917 


.93083 


.06763 


.00153 




PHOTOS EXP 


.93087 


.06913 


.93087 


.06747 


.00166 





.0046 





.0041 


.0255 



exact TAUOLA is compared with single-photon mode of PHOTOS. One can notice that the 
difference between the single-photon emission mode in PHOTOS and TAUOLA is similar to 
the difference between exponentiated mode of PHOTOS and TAUOLA. This indicates that the 
dominant source of differences is due to non-leading 0{a) terms that are missing in PHOTOS 
(and present in full 0(a) ME of TAUOLA). Predictions of the single-photon mode of TAUOLA 
are more precise than PHOTOS predictions, even with exponentiation. Nevertheless, it must be 
realized that as TAUOLA does not implement exponentiation, it is limited to the generation of 
(at most) single hard photon emission. For cases where the topologies of configurations with 
two or more photons are important, the exponentiated algorithm implemented by PHOTOS will 
be more appropriate. 

5.3 Tests of PHOTOS at very high energies 

Because of numerical instabilities, the use of PHOTOS at LHC energies of 14 TeV was techni- 

cally limited before, by a restriction on minimal ratio (M being the mass of the decaying 
particle, E™" the minimal energy of generated photons). The phase space for generation had to 
be significantly reduced and photons of even moderate energies could not be generated. 

We verified that the exponentiated version of the algorithm, armed with the new kinematics- 
correction routine implemented in PHOTOS 2.13, allows it to be used to generate bremsstrahlung 
in decays of very energetic particles at the energy scales that indeed require very low cut on the 
minimal photon energy. 
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Figure 7: Comparison of full 0(a) in TAUOLA and exponentiated version of PHOTOS for 
x~ — i> /j^VT;V^(y) channel, iStest = 0.05 GeV, testl is used. The plot below presents the distribu- 
tion of invariant mass of Vi;V^Yfrom TAUOLA (in red, or darker-grey) and PHOTOS (in green, 
or lighter- grey); the black curve is the ratio of the two normalized distributions. 



GENERATOR 
n photons 


Branchi 



ng ratio 
1 


Max SDP 
1 1 


TAUOLA 


0.98916 


0.01084 




PHOTOS EXP 


0.98935 


0.01065 





0.0047 




18 



We performed studies using KKMC as a "host" generator for PHOTOS to produce leptonic 
Z-decays where the mass of the produced Z/y intermediate state was of order 2 TeV. The value 

^min 

of the infrared cut-off parameter was lowered down to the value of ~ 10^^, which allowed 
the generation of relatively soft photons, and high-statistics runs were completed without en- 
countering numerical stability problems. Again, we have found excellent agreement between 
the exponentiated version of PHOTOS and the exponentiated 0{a^) KKMC. As usual, com- 
plete numerical results can be found on the web page llISI . 

5.4 Recent developments : semileptonic K decays and universal interfer- 
ence weight 

Recently our attention was drawn to a use of PHOTOS in an estimation of the effects of QED 
radiative corrections in high-precision measurements of K, B and D meson decays. The system- 
atic error of results given by PHOTOS is important for measurements of elements of Cabibbo- 
Kobayashi-Maskawa matrix of quark mixing, because QED corrections affect the acceptance. 
The insufficient precision of PHOTOS observed in ref. PUIl was of no surprise: for the case of 
semileptonic K decays, PHOTOS was prepared for use in the "crude" mode only, suitable for 
full phase space coverage of bremsstrahlung for detector studies, but not for shape estimations. 
The most important missing terms in calculations performed by PHOTOS in this case were 
those from QED interference. 

Thanks to the technical development performed in PHOTOS, and partially described already 
in this paper, it was rather straightforward to implement the universal interference weight for 
any decay channel, as expressed by formula (17) of Q: 



Wr 



multi 



qi-k 



+ 2 



2 

2 - " 



qi-k 



(1) 



+ .. 



where Qi, Qi,... denote the charges of the decay products, q\, qi,... denote the momenta of 
the decay products, k denotes the energy of photon, and a summation is performed over photon 
polarization states denoted by £. 

Earlier, the interference weight in PHOTOS was calculated from internal angular variables 
and not from four- vectors; therefore such universal implementation was not trivial. The impact 
of the new interference weight on the photon angle distributions in ^ ^ A'V7c(y) process is 
presented in Fig.lU It seems that the majority of discrepancies between PHOTOS and the first- 
order ME generator KLOR have been removed in this way^^. The complete study will however 
be documented elsewhere iHDl . 
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For the complete comparison, see Fig. 8 of BOI and our comparison web page 11381 . 
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Figure 8: Impact of the universal interference weight in PHOTOS on predictions for the angular 
distribution in ^ ^ A'V(y) events. The distribution of the cosine of the angle between the charged 
lepton and the photon generated by older version of PHOTOS, without the interference weight 
(plotted as black dots), and PHOTOS 2.13 with universal interference weight (plotted as red 
line). This figure should be compared with figure 8b from [40J. The majority of discrepancies 
between PHOTOS and 0(a) generator KLOR seems to be removed by using the universal 
interference weight. 
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6 Summary and Outlook 



We have presented numerical results from a recently upgraded version of the PHOTOS Monte 
Carlo algorithm for radiative corrections in decays of resonances. In the case of W and Z 
decays, we found the precision of the simulations in complete multiple-photon final states at 
the level of (or better than) 0.1%. This conclusion originates from numerous comparisons 
with precision Monte Carlo programs: KORALZ, KKMC and WINHAC for simulation of 
radiative corrections in Z and W decays. The reliability of those programs - from theoretical 
and technical sides - was well established by numerous applications at LEP time, documented 
in a multitude of publications on comparisons with the experimental data (in the case of Z 
decay), as well as other programs and calculations. 

The results we presented prove numerically that PHOTOS may be used with confidence to 
calculate bremsstrahlung in Z and W leptonic decays with very high precision. In the case of 
Higgs boson decays and leptonic decays of x's, this precision tag is valid as well, even if the 
proof is not as solid, since it relies on the test at first order only. 

This conclusion cannot be extended to other decay channels, without further work involving 
a study of process-dependent matrix elements. Nonetheless, precision in the general case is also 
better now. 

Finally, let us stress that the mathematical side of the multiple-photon algorithm presented 
here is missing a rigorous proof. Present-day applications do not justify such an effort. Wher- 
ever it was possible, we have verified PHOTOS correctness and found that the performance was 
excellent. In other cases we were limited by availability of exact matrix elements for single- 
photon emission. 

A more thorough explanation of the mathematical side of the algorithm may, however, be 
interesting in the future, especially having in mind applications in QCD. Note that algorithms 
of PHOTOS do not rely on the conformal symmetry of soft photon factors and the phase space 
of masless particles, which is crucial in the construction of exclusive-exponentiation algorithms 
as used by KKMC, WINHAC or KORALZ. 

From its origin, the algorithm of PHOTOS was developed as FORTRAN77 code, even 
though a C-i-i- implementation was already completed several years ago. Since there was 
no standard for C-I-+ event record at that time, it was not published. Recently, the experimental 
communities expressed interest in using PHOTOS with the HepMC event record P31l . This 
certainly is a valid option for future development. 
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A New features of PHOTOS version 2.13 



In this appendix let us briefly summarize options and switches available in version 2.13 of PHO- 
TOS, so it can serve as a reference guide for librarians and end-users. The complete, technical 
description of the algorithm, with its recent improvements, will be presented elsewhere lfT2lfT3l . 

By default, the new options are not active and PHOTOS 2.13 behaves much like PHOTOS 
2.0 from 1994 |l2|. It also does not require any changes in old interfaces. In most applications it 
will produce results that differ little with respect to previous versions, with the exception of the 
case of interference effects which are now simulated. As in previous versions, the initialization 
parameters, which are subject to user modifications, should be configured in the subroutine 
PHOCIN. 

A.l Standard settings 

In Table |3] we summarize all switches used in standard PHOTOS setups. The names of the 
options that are new (or need special care) in version 2.13 are typed in bold font. 

Exponentiated multiple-photon mode 

The exponentiated mode has been implemented in version 2.09 of PHOTOS; it may be activated 
using the I EXP flag. By eliminating the problem of negative probabilities we were able to lower 
the default value for XPHCUT parameter to 10^^ for that mode. Nevertheless, multiple iteration 
of emission kernel puts higher demands on numerical precision: aggregating rounding errors 
and initially insignificant numerical deficiencies in energy-momentum conservation of the input 
data, may lead to significant problems with numerical stability. Therefore a new option for 
kinematics corrections of the data in the event record was implemented^^. The new option, 
activated with CALL PHC0RK(5) is active by default for the exponentiated mode. It corrects 
momenta four-vectors of all particles taking part in photon generation: at first, the energies of 
all daughter particles are corrected in such a way that eJ = pj +Mf; then the four-momentum 
of the mother-particle is recalculated to ensure energy-momentum conservation. 

Universal interference weight 

Before version 2.13 the QED interference was implemented in PHOTOS for dedicated decay 
channels only, namely neutral particle decaying into a pair of opposite-charge particles For 
other channels the interference effects were not taken into account so far. 

Starting from version 2.13 a new, universal algorithm calculating interference effects for 
all processes is present in PHOTOS. By default, as in previous versions, this algorithm (im- 
plemented internally as the Monte Carlo interference weight) is active and may be controlled 

'^This is the extension of the kinematic correction schemes described in subsection 7.2.1 of H LI ■ 

'"*This requirement was relaxed slightly in PHOTOS 2.07. Since then the pair of decay products could be of 

different masses. See subsection 7.2.2 of 1311 for details. 

'^Except for leptonic W decays, where a dedicated interference weight was implemented in version 2.07 (see 

subsection 7.2.4 of ET\ ). 
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Table 3: Standard PHOTOS switches and configuration parameters, (exp) denotes the switches 
and values that are specific to exponentiated mode. 



Subroutine 


Parameter 


Default 


Description 


Options for photon multiplicity selection 


PHOCIN 


ISEC 


.TRUE. 


Switch to enable double emission 


PHOCIN 


ITRE 


.FALSE. 


Switch to enable triple/quartic emission 


PHOCIN 


lEXP 


.FALSE. 


Switch to enable exponentiated, multiphoton mode 


Options related to interference 


PHOCIN 


INTERF 


.TRUE. 


Switch to enable universal interference weight 


PHOCIN 


IFW 


.TRUE. 


Switch to enable dedicated interference weight 
for W lv{y) decays 


PHOCIN 


PINT 


2.D0 


Maximum interference weight. 
The value may need to be adjusted to 2"^^ , where 
n is the maximal multiplicity of charged decay 
products of any elementaiy branching 


Other tec 


mical options 


PHOCIN 


CALL PHCORK(n) 


n=l 
n=5 (exp) 


Switch to select kinematics correction mode: 
n=l: no correction 
n=2: corrects energy from mass 
n=3: con^ects mass from energy 
n=4: coiTccts energy from mass for particles 
up to 0.4 GeV mass, for heavier ones corrects mass 
n=5 : most complete correction for mother 
and daughter particles; required for (exp) 


PHOCIN 


XPHCUT 


0.01 
lD-7 (exp) 


Infrared cut-off parameter: minimal photon energy 
expressed as fraction of decaying particle mass 


PHOCIN 


IFTOP 


.TRUE. 


Switch to enable emission in the 
hard process gg{qq) — > tt 
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Table 4: Advanced PHOTOS switches and configuration parameters 



Subroutine 


Parameter 


Default 


Description 


PHTYPE 


IFOUR 


.TRUE.AND.(ITRE) 


Switch to enable quartic emission 


PHOCIN 


EXPEPS 


lD-4 


Residual probability for multiple iteration, 
"stop criteria" for crude generation 



Selection of weighting algorithm variant 



PHOCIN 
PHOMAK 
PHOPRE 
PHOENE 
PHOINT 






Select the variant of weighting algorithm 
code. May be altered by proper commenting 
of code blocks denoted with VARIANT A 
and VARIANT B comments. 
For use by experts only. 



using INTERF flag. To use this new interference weight, it is required that the value of the FINT 
parameter be adjusted in such a way that the maximum weight in the Monte Carlo rejection is 
increased to 2"^^, where n denotes the maximal expected multiplicity of charged particles in 
decays processed by PHOTOS. 

An omission in adjusting this parameter often results in PHOTOS stopping with an error 
message indicating too high a value of summary weight. On the other hand, putting the value 
of FINT too high may significantly increase the CPU time consumed by PHOTOS because 
of Monte Carlo rejection not being effective. An option for having the actual value of FINT 
automatically calculated for each event is being considered for implementation in the future. 

A.2 Advanced settings 

In Table m we summarize the switches and options dedicated for advanced users. We strongly 
recommend to contact the authors before these settings are altered. 

At an intermediate step of development, a more effective Monte Carlo weight for crude 
probability in interference calculation was implemented in PHOTOS 2.07: instead of doubling 
the crude distribution for photon emission from each charge, a flat, parallel channel was added. 
Although the results were satisfactory, this solution turned out to suffer from problems with 
stability when used with universal interference weight for more that three-body decays. We 
therefore reverted the code to the original implementation. For the purpose of preserving this 
interesting technical development and maintaining the compatibility with PHOTOS versions 
2.07-2. 12, released in 2003-2005, we leave this optimized code as an option to be activated by 
advanced users only. The two options are marked in the code by VARIANT A and VARIANT B 
comments, the former being currently active, the latter being commented-out. Due to elimina- 
tion of numerical instabilities, we were able to lower the value of XPHCUT for variant B to 10^^ 
and perform important tests. Such tests were not possible in the default variant A, because the 
value of XPHCUT has to be kept at 10"^. 
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